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Uniform electron gases: III. Low-density gases on three-dimensional spheres 

Davids Agboola,^ Anneke L. Knol,^ Peter M. W. Gill/’[^and Pierre-Fran^ois Loos^ p^ 

Research School of Chemistry, Australian National University, ACT 2601, Australia 

By combining variational Monte Carlo (VMC) and complete-basis-set limit Hartree-Fock (HF) calculations, 
we have obtained near-exact correlation energies for low-density same-spin electrons on a three-dimensional 
sphere (3-sphere), i.e. the surface of a four-dimensional ball. In the VMC calculations, we compare the 
efficacies of two types of one-electron basis functions for these strongly correlated systems, and analyze the 
energy convergence with respect to the quality of the Jastrow factor. The HF calculations employ spherical 
Gaussian functions (SGFs) which are the curved-space analogs of cartesian Gaussian functions. At low 
densities, the electrons become relatively localized into Wigner crystals, and the natural SGF centers are 
found by solving the Thomson problem (i.e. the minimum-energy arrangement of n point charges) on the 
3-sphere for various values of n. We have found 11 special values of n whose Thomson sites are equivalent. 
Three of these are the vertices of four-dimensional Platonic solids — the hyper-tetrahedron {n = 5), the 
hyper-octahedron [n = 8) and the 24-cell (n = 24) — and a fourth is a highly symmetric structure (n = 13) 
which has not previously been reported. By calculating the harmonic frequencies of the electrons around 
their equilibrium positions, we also find the first-order vibrational corrections to the Thomson energy. 


I. INTRODUCTION 


In a recent series of papers we have shown that 
the behavior of electrons in the (flat) Euclidean space 
is surprisingly similar to the behavior in the (curved) 
Zl-dimensional manifold S^, the surface of a (U -F 1)- 
dimensional ball. ^ By e xploiting this similarity betw een 
electrons on a lin^iSin] and electrons on a rin^SHHHl^ 
have constructed a new type of correlation functional 
(the generalized local-density approximation) for density 
functional theory (DFT) calculations and we have shown 
that this new functional yields accurate corr elat ion ener¬ 
gies in a variety of one-dimensional systems.^^ 

However, our ultimate goal is to construct improved 
functionalsiSl for three-dimensional (3D) systems and, to 
this end, we seek accurate correlation energies from the 
spherical space (henceforth called a glome), which is 
the surface of a four-dimensional (4D) ball. Electrons 
in or enjoy three degrees of translational freedom 
but the properties of n electrons on a glorne ( n-glomium) 
have hitherto received little attention.EE3Hl2l 

An n-glomium is defined by the number n of electrons 
and the glome radius R. Its electron density is 

p = n/{2TT'^R^), ( 1 ) 

but it is often measured via the Wigner-Seitz radiu^^ 


2n J 


R, 


( 2 ) 


which measures the average distance between neighbour¬ 
ing electrons. High-density systems (which are weakly 
correlated) have small values while low-density sys¬ 
tems (which are strongly correlated) have large val¬ 
ues. In the present study, we focus our attention on low- 
density glomiums with 2 < n < 48. 
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In the low-density regime, the Coulomb energy (which 
decays as R~^) dominates over the kinetic energy (which 
decays as R~^) and the n electrons localize onto partic¬ 
ular points on the glome that minimize their (classical) 
Coulomb repulsion. These minimum-energy configura¬ 
tions are called Wigner crystalJi^ and, if all of its sites 
are topologically equivalent, we will call it a uniform lat¬ 
tice. 

The work is organized as follows. In Sec. |llj we study 
the Thomson problem on a glome and discuss the uni¬ 
form solutions. In Sec. m we calculate the harmonic vi¬ 
brational energy of the electrons as they oscillate around 
their lattice positions. In Secs. |l^ and respectively, 
we report Hartree-Eock, near-exact and correlation en¬ 
ergies of n-glomium at various densities. Unless other¬ 
wise stated, all energies are reduced (i.e. per electron). 
Atomic units are used throughout. 


II. THE THOMSON PROBLEM 

What arrangement of n unit point charges on a unit D- 
sphere minimizes their classical Coulomb energy? This 
generalizes a question posed by J. J. Thomson as he de¬ 
vised the “plum pudding” model of atomic structure.!^ 
Although the model was abandoned long ago, the Thom¬ 
son problem continues to intrigue mathematicians and 
has resurfaced in many fields of science: surface or¬ 
dering of liquid metal drops confined in Paul traps,!^ 
fullerenes-like molecules,^ arrangements of protein sub¬ 
mits on sphe rical virusej^ or multielectron bubbles in 
liquid helium.l^^ 

The Thomson problem on a 1-sphere (i.e. a ring) is 
trivial and the solutions consist of charges uniformly 
spaced around the ring. The problem on a 2-sphere is 
challenging and, although it has been studied numeri¬ 
cally up to large values of mathematically rig¬ 

orous solutions^S have been established only for n € 
{2,3,4,5,6,12}. 

The Thomson problem on a 3-sphere (i.e. a glome) 
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TABLE 1. Non-reduced Thomson energy eg and principal moments of inertia Ik of the Thomson lattices on a unit glome. 


_n _£o_ 

2 1/2 

3 3/^ 

4 6/^873 

5 

6 9/V^ + 6/V3 

8 24/^+ 4/2 

10 ^ 10 
2sin(7r/5) ^ 2sin(27r/5) 

ifl 6 12 12 24 ^ 12 

2c yss V2 + 2s2 V2 - s2 V4- s2 

13 __ - - _ 

.id,4 2 cos(A:7r/13) cos(5fc7r/13) 

24 96/1 + 72/d2 + 96/d3 +12/2 

24 240 48 96 96 

48 TT + ^ + / + • + ' ■ , 

2 V2 V2±V2 ^2 + 72/2 ^2 + 76/2 

^ 96 96 

^2±(d3-l)/2 y'2±(d3 + l)/2 


h 

h 

h 

h 

0 

2 

2 

2 

3/2 

3/2 

3 

3 

8/3 

8/3 

8/3 

4 

15/4 

15/4 

15/4 

15/4 

9/2 

9/2 

9/2 

9/2 

6 

6 

6 

6 

15/2 

15/2 

15/2 

15/2 

6(1 Ac'*) 

6(1 Ac"*) 

6(1 A s") 

6(1 A s^) 

39/4 

39/4 

39/4 

39/4 

18 

18 

18 

18 

36 

36 

36 

36 


^ c = cos s = sin 0 and 0 = 0.7935536685... 


seeks the global minimum eg = l^(xo) of 

n 

i<3 

where x describes the positions of the n charges on the 
glome and is the Euclidean distance between charges i 
and j, measured th rough the unit glome. It has attracted 
much less attentioiP^El] than the D = 2 problem. 

We performed a numerical study to determine the val¬ 
ues of n < 50 for which the Thomson minimum-energy 
configuration on a glome is uniform. Although it would 
be more suitable to use a global optimization, because we 
consider relatively small numbers of electrons, we adopt 
the following computational strategy: for each n, we gen¬ 
erated randomly at least 1000 distinct initial structures 
and minimized their energy using local optimization al- 
gorithms, as implemented in the Mathematica software 
package.l^Our numerical experiments indicate that there 
are eleven uniform lattices. Their energy Sq and principal 
moments of inertia are listed in Table Q] and their carte¬ 
sian coordinates are in supplementary material.l^ The set 
{2,3,4, 6,8,12, 24} of n values for uniform D = 2 Thom¬ 
son latticeJi^ is a subset of the D = 3 set but we do not 
understand this. 

The principal moments of inertia of a lattice indicate 
the degree of its symmetry. Generalizing the standard 
notation for 3D structureswe define a “hyperspherical 
top” as a 4D structure with four equal moments; the 
extremely symmetrical n = 5, 6, 8, 10, 13, 24 and 48 
lattices are of this type. We define a “spherical top” as 
a lattice in which three of the four moments are equal; 


the n = 2 (prolate) and n = 4 (oblate) lattices are of 
this type. We define a “symmetric top” as a lattice in 
which the moments form two pairs; the n = 3 and n = 12 
lattices are of this type. 

The glome lattices for n = 2 (a diameter) and n = 3 
(equilateral triangle) are the same as on a 1-sphere and 2- 
sphere. The glome lattice for n = 4 (regular tetrahedron) 
is the same as on a 2-sphere. The n = 5 lattice is a reg - 
ular hyper-tetrahedron (also called a regular simplejP^, 
a 4D Platonic solid with ten equal side lengths. The 
n = 6 lattice is the union of an equilateral triangle in 
the wx-plane and another such triangle in the yz-plane. 
The n = 8 lattice is a hyper-octahedron (or 16-cell), a 
4D Platonic solid with vertices at ±1 on each of the four 
cartesian axes. The n = 10 lattice is the union of a regu¬ 
lar pentagon in the r(;a;-plane and another such pentagon 
in the yz-plane, while the n = 12 lattice is the union of 
two perpendicular triangular prisms. The n = 13 lattice 
is peculiar to the D — 3 Thomson problem and, to the 
best of our knowledge, has not been previously described. 
The n = 24 lattice is the 24-cell, a 4D Platonic solid with 
no analogue in 3D. It is the union of a hyper-octahedron 
and a hyper-cube. The n = 48 lattice is peculiar to the 
D = 3 Thomson problem. 


III. HARMONIC VIBRATIONAL ENERGY 

The energy of an n-electron Wigner crystal on a 
glome can be estimated by solving the Schrodinger equa- 
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TABLE II. Reduced Thomson energy Eo, harmonic vibra¬ 
tional energy Ei, near-exact energy E, Hartree-Fock energy 
Ehf and correlation energy E^ (in mi5h) for various n and Vs 




Eg 

20 

50 

100 

150 



Eo 

16.633 

6.653 

3.327 

2.218 



Eq + El 

23.068 

8.281 

3.902 

2.531 

n = 

■■ 2 

E 

23.928 

8.385 

3.924 

2.540 



Ehf 

24.911 

8.795 

4.100 

2.643 



-E, 

0.983 

0.410 

0.176 

0.103 



Eq 

33.557 

13.423 

6.711 

4.474 



Eq + El 

41.074 

15.324 

7.384 

4.840 

n = 

: 3 

E 

41.783 

15.405 

7.400 

4.847 



Ehf 

43.811 

16.079 

7.664 

5.000 



-Be 

2.028 

0.674 

0.264 

0.153 



Eq 

48.507 

19.403 

9.701 

6.468 



Eq + El 

57.190 

21.600 

10.478 

6.890 

n — 

: 4 

E 

57.155 

21.550 

10.461 

6.882 



Ehf 

59.886 

22.358 

10.757 

7.046 



-Be 

2.731 

0.808 

0.296 

0.164 



Bo 

62.009 

24.804 

12.402 

8.268 



Eq + El 

71.916 

27.310 

13.288 

8.750 

n — 

: 5 

E 

71.038 

27.119 

13.238 

8.728 



Ehf 

74.240 

27.985 

13.537 

8.888 



-Be 

3.202 

0.866 

0.299 

0.160 



Bo 

75.564 

30.226 

15.113 

10.075 



Eq + El 

85.823 

32.821 

16.030 

10.575 

n — 

: 6 

E 

85.406 

32.682 

15.986 

10.553 



Ehf 

88.283 

33.530 

16.293 

10.720 



-Be 

2.877 

0.848 

0.307 

0.167 



Eq 

99.390 

39.756 

19.878 

13.252 



Eq + El 

110.951 

42.681 

20.912 

13.815 

n — 

: 8 

E 

110.614 

42.620 

20.915 

13.823 



Ehf 

112.893 

43.234 

21.115 

13.927 



-Be 

2.279 

0.614 

0.200 

0.104 



Eq 

122.336 

48.934 

24.467 

16.311 



Eq + El 

133.974 

51.879 

25.508 

16.878 

n — 

10 

E 

133.522 

51.767 

25.485 

16.874 



Ehf 

136.154 

52.479 

25.727 

17.002 



-Be 

2.632 

0.712 

0.242 

0.128 



Bo 

143.339 

57.336 

28.668 

19.112 



Eq + El 

155.615 

60.441 

29.766 

19.710 

n — 

12 

E 

154.713 

60.269 

29.740 

19.707 



Ehf 

157.598 

61.001 

29.970 

19.822 



-Be 

2.885 

0.732 

0.230 

0.115 



Bo 

153.600 

61.440 

30.720 

20.480 



Eq + El 

165.909 

64.554 

31.821 

21.079 

n — 

13 

E 

164.804 

64.322 

31.767 

21.061 



Ehf 

167.947 

65.126 

32.033 

21.196 



-Be 

3.143 

0.804 

0.266 

0.135 



Eq 

252.272 

100.909 

50.454 

33.636 



Eq + El 

265.600 

104.280 

51.647 

34.285 

n — 

24 

E 

264.917 

104.138 

51.624 

34.280 



Ehf 

267.464 

104.793 

51.831 

34.387 



-Be 

2.547 

0.655 

0.207 

0.107 



Eq 

425.792 

170.317 

85.158 

56.772 



Eq + El 

439.690 

173.833 

86.401 

57.449 

n — 

48 

E 

438.667 

173.681 

86.406 

57.474 



Ehf 

441.542 

174.334 

86.581 

57.547 



-Be 

2.875 

0.653 

0.175 

0.073 

n — 

oo 

-Be 

6.839 

3.607 

2.080 

1.499 


tion in the harmonic potential 

V' 2 (x) = y(xo) -h ^(x-Xo)"^ • H • (x-Xo), (4) 


where H is the 3nx3n second-derivative (Hessian) matrix 


_ a^H(x) 

dtidtj 


(5) 


and the ti are suitable tangential coordinates. The square 
roots of the Hessian eigenvalues are the harmonic fre¬ 
quencies uji and one can then write 


^w=§ + ^fo^^+0(R->) (6) 

= Eo + Ei+0{R-^). (7) 

One finds that, for n > 3 particles in S^, exactly 3n — 6 
of the Hessian eigenvalues are non-zero and six vanish 
because they correspond to rotations on the glome. This 
is analogous to the familiar 3n — 6 rul^^ for non-linear 
molecules vibrating in K^. Numerical values of Eq + Ei 
for a range of n and are presented in Table |TTj 


IV. HARTREE-FOCK ENERGIES 

We now turn to the ab initio treatment of n spin-up 
electrons on a glome, i.e. ferromagnetic n-glomium. We 
have performed Hartree-Fock (HF) calculationJ^in a ba¬ 
sis of s-type spherical Gaussian functions (SGFs)P^ 

where A G is the center of the SGF, a is the exponent 
and Ii is a modified Bessel function.ISS An SGF behaves 
like a Gaussian near A and is therefore a natural basis 
function for describing a localized electron. Moreover, 
the product of two SGFs is a third SGF which make 
them computationally attractive.!^ All the required one- 
and two-electron integrals can be found in Ref. IT71 

Adoptin g th e “Gaussian lobe” philosophy introduced 
by Whitten^ many years ago, we use off-center s-type 
SGFs to mimic SGFs of higher angular momentum. The 
basis set consists of a grid of s-type SGFs with same expo¬ 
nent a clustered around each Thomson site (see Fig.[^. 
The complete basis set (GBS) HF energy is obtained by 
successively adding Level 0 (LO), Level 1 (LI), Level 2 
(L2) and Level 3 (L3) functions to the basis set. In each 
calculation, we optimize the SGF exponent a and the 
LO/Ll distance 8 using the Newton-Raphson optimiza¬ 
tion procedure. Our target accuracy was 1 microhartree 
(/j,Ah) per electron. The resulting HF energies for a range 
of n and are shown in Table [TTl 

Table [TTT| reports HF energies of n-glomium as the basis 
set is gradually improved. SGFs are optimal for localized 
electronic systems but become less efficient as the density 
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TABLE III. HF energies (in mEh) for ferromagnetic n-glomium with Vs = 20. There are M basis functions per Thomson site. 


Basis set 





Number of electrons n 





Level 

M 

2 

3 

4 

5 

6 

8 

10 

12 

13 

24 

48 

0 

1 

24.983 

43.939 

60.016 

74.277 

88.345 

112.911 

136.249 

157.626 

167.991 

267.468 

441.543 

1 

7 

24.917 

43.816 

59.890 

74.253 

88.290 

112.899 

136.169 

157.606 

167.951 

267.464 

441.542 

2 

19 

24.911 

43.811 

59.886 

74.244 

88.286 


136.156 

157.605 

167.947 

5? 

?! 

3 

27 

5) 







5) 


?? 

!? 


<t> 



FIG. 1. Three-dimensional SGF grid used at each Thomson 
site in HF calculations. The Level 0 (LO), Level 1 (LI), Level 
2 (L2) and Level 3 (L3) functions are represented in purple, 
red, green and blue, respectively. S is the shortest distance 
(around the glome) between two grid points. 

increases so the convergence for Cg > 20 is always at least 
as fast as for Vg = 20. We therefore show results for 
Tg = 20, the most challenging case. 

For a given value of Cg, the minimal-basis (LO) expo¬ 
nent a grows, i.e. the electrons become more localized, 
as n increases. The results of Table m show that L2 
achieves /iifh accuracy for all n values and, indeed, LI 
suffices for the largest n values. It is well known that, on 
a 2-sphere, the number of nearest neighbors around an 
electron approaches six (hexagonal lattice) for large 
Similarly, on a glome, the number of nearest neig hbors 
approaches eight (body-centered cubic lattice).^ Thus, 
for large n, the density around each electron becomes 
approximately isotropic, and the L2 and L3 functions 
become largely superfluous. 

V. NEAR-EXACT ENERGIES 

To obtain near-exact energies for electrons on a glome, 
we have performed variational Monte Carlo (VMC) 


calculations.!^ The trial wave function is of the form 

T'T = «'oe'^, (9) 

where 410 is a Slater de termi nant of either SGF^i^ or 
hyperspherical harmonic^^^!^ (HSHs) 

Yumix^d.cj)) = C'f+^(cosx)sin^xWm(6»,(^). (10) 

is a Gegenbauer polynomial and is a spherical 
harmonic.!^ The Jastrow factor J is a symmetric function 
of the interelectronic distances containing two- body 
(2B), three-body (3B) and four-body MB) terms.l^^I^ 
More details will be reported elsewhere.S3 

At low densities, the energy minimization procedure is 
unstable and the parameters of the Jastrow factor were 
therefore optimized by variance minimization using New¬ 
ton’s method.!^2!!^In all calcul ations , the statistical error 
obtained by reblocking analysi^^^ESl jg under 1 /rifh- For 
small numbers of electrons, comparisons with extrapo¬ 
lated full configuration interaction (FCI) calculation^Smil 
indicate that our VMC energies have sub-^Ah accuracy. 
They are reported in Table|ll]for various rg and n values. 

Because we have observed that many-body effects are 
more important for small rg, we have studied the con¬ 
vergence of the energy for rg = 20 and various n values 
as 2B, 3B and 4B terms are successively included. We 
found that the inclusion of 3B terms systematically low¬ 
ers the energies by up to 20 /rAh per electron but that the 
inclusion of 4B terms offers less than 1 /iAh per electron. 
We therefore eschewed 4B terms in the calculations with 

rg > 20. 

When SGFs were used, the determinant corre¬ 
sponds to a HF Level 0 calculation (i.e. a single SGF 
on each Thomson site) but the value of the SGF expo¬ 
nent was optimized to minimize the VMG energy. We 
found that the SGF basis is superior to the HSH basis 
for n = 8, 10, 12 and 24, while the two basis sets yield 
identical energies for the other cases. 8-, 10-, 12- and 
24-glomium are “open-shell” systems, i.e. the highest oc¬ 
cupied shell of HSHs is partially occupied and there are 
several low-lying determinants with signihcant weights. 
In contrast, SGFs at the Thomson lattice sites naturally 
describe the localized electrons and are particularly well 
suited to these systems. For 48-glomium, computational 
limitations precluded full exponent optimization and the 
results in Table m were therefore obtained with HSHs. 
They are probably less accurate than the other energies. 
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Ec (niEh) 



FIG. 2. Reduced correlation energy Ec (in miJh) as a Inac¬ 
tion of n for varions Vs- 


VI. DISCUSSION 

By taking the difference between the CBS-HF ener¬ 
gies of Sec. |I^ and the VMC energies of Sec. |Vj we 
have obtained the near-exact correlation energies Ef. of 
n-glomium for various and n values. Our results are 
reported in Table [l^ 

While we have shown that, in ID, the correlation en¬ 
ergy is a smooth and monotonic function of the 
situation is rather different in 3D. As shown in Figure 

(where have plotted Ec as a function of n), the re¬ 
duced correlation energies do not change monotonically 
as n increases. Instead, they initially increase and reach 
a maximum at n = 5 or n = 6. Beyond n = 5, they os¬ 
cillate and tend to decrease slowly with n, especially at 
very low densities. The oscillations are probably due to 
“shell effects” which originate from partially-filled energy 
level in open-shell systems (see above). Such shell effects 
are also observed in 2D For fixed n, our numerical re- 

_Q /2 

suits show that E^ decreases as rs for large r^. This 
is expected due to the cancellation of the leading term 
(proportional to in the exact and HF energies ex¬ 
pansion at large rg (see Eg. ([^). We have also reported 
the correlation energies^^l*^ of the jellium model (which 
corresponds to n — > oo) in Table [TT} Exact jellium ener¬ 
gies are available i n man y previous papers^SHaSl but fully 
relaxed HF energie^^^*^ are rare. As a result, most of the 
jellium correlation energies in the literature are based on 
unrelaxed HF energies and are consequently significantly 
larger than ours. 

The harmonically corrected Thomson energy (Eq+Ei) 
is usually higher than the exact energy but n = 2 and n = 
3 are exceptional cases. At very low densities, however, 
it always approximates the exact energy well. Including 
the first enharmonic correc tion E 2 would probably yield 

even better estimates 


VII. CONCLUSION 

The goal of this work was to generate benchmark cor¬ 
relation energies for the development of improved corre¬ 
lation functionals for DFT calculations on 3D systems. 


To achieve this, we have studied the correlation energies 
of low-density spin-polarized electron gases on a glome. 

First, we looked at the Thomson problem on a glome, 
using numerical optimization algorithms to locate uni¬ 
form lattices, i.e. those in which all the lattice sites are 
equivalent, for n < 50. We found eleven uniform lattices. 
We also noted that the structure of three of these uniform 
lattices correspond to well-known 4D Platonic solids: the 
hyper-tetrahedron (n = 5), the hyper-octahedron {n = 8) 
and the 24-cell (n = 24). Moreover, we have pointed out 
the highly symmetric case of n = 13 and we stressed 
that this polychoron has not been previously described 
anywhere in the literature. By taking into account the 
quantum oscillation of the electrons around their equi¬ 
librium positions, we obtained the harmonic vibrational 
contribution to the (classical) Thomson energy. As ex¬ 
pected, the sum of the Thomson and harmonic energies 
is a very good approximation of the exact energy at very 
low density. 

Moreover, by systematically increasing the number of 
s-type SGF basis functions around each Thomson site, we 
obtained the CBS HF energy of n-glomium for a range 
of densities (20 < rg < 150). In general, the convergence 
analysis reveals that only L2 calculations are required 
to converge the HF energies to microhartree accuracy. 
However, we note that as we move into the high-density 
regime, more s-type SGF functions will be needed to 
reach the CBS limit. In this regime, HSHs might con¬ 
stitute a more suitable one-electron basis set. We will 
investigate this in a forthcoming paper. 

The near-exact energies were obtained using highly- 
accurate VMC stochastic calculations using HSH and 
SGF one-electron basis sets. We have shown that 4B 
terms in the Jastrow factor have insignificant effect on 
the energy. The energies obtained using both basis sets 
are in very good agreement, except for 8-, 10-, 12- and 
24-glomium where the SGF basis is superior. 

The present work is a significant step towards the 
construction of correlation functionals for molecules and 
solids. The next step is to generate accurate correlation 
energies within the high density regime, and for partially- 
polarized systems. This work is underway and our results 
will be reported elsewhere. 
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